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Abstract 



How fast does the quark condensate in QCD-like theories vary as a function of Nf 
is inferred from real QCD using chiral perturbation theory at order one-loop. A sum 
rule is derived for the single relevant chiral coupling-constant, Lq. A model independent 
lower bound is obtained. The spectral function satisfies a Weinberg-type super convergence 

, relation. It is discussed how this, together with chiral constraints allows a solid evaluation 
of Lq, based on experimental tttt — KK S-wave T-matrix input. The resulting value of 
Lq is compatible with a strong Nf dependence possibly suggestive of the proximity of a 

. chiral phase transition. 

as : 
o . 

Q\ • 1 Introduction 
OV 

By analogy with recent results obtained in supersymmetric theories [Q, one expects that 
Qh| QCD-like theories will undergo a number of phase transitions at zero temperature upon 
varying Nf, the number of different flavour fermions, at fixed number of colours (N c = 3 in 
the following). If the number of fermions is large Nf > \-N c the theory has no asymptotic 
freedom and no confinement. Decreasing Nf below -^N c one encounters a conformal 
phase as indicated by the fact that the j3- function at two loops has a zero[[|. Assuming 
^ . the fermions to be all massless, the chiral SU (Nf) x SU(Nf) symmetry of the QCD action 
remains unbroken in this phase. If one further decreases Nf to small values, Nf = 2 — 3, 
then QCD is in a confining phase in which the chiral group is spontaneously broken to 
SU(Nf). It is generally believed that, in this phase, the quark condensate is non- vanishing 
and largeQ. At larger Nf, there could exist different phases of chiral symmetry breaking. 
A physical picture of such phases is proposed in ref.0. An interesting open question 
concerns the value of Nj nt for which a chiral phase transition takes place. Recent lattice 
results suggest that a transition could occur for Nf as small as four|5|. In an instanton 
vacuum model, the quark condensate ceases to be non-vanishing for Nf of the order of 
five|J, while another theoretical model obtains a much larger value [j?]]. 

In this paper, we use the fact that nature solves ordinary QCD in order to extract an 
information on how the quark condensate varies with Nf for small values of Nf. More 

1 Address after Sept. 1: MIT, Center for Theoretical Physics, 77 Massachussets Av., Cambridge MA 
02139-4307 

2 Experimental verification of this conjecture necessitates specific and very precise data. This question 
is discussed in ref. 131. 
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specifically, if it can be shown that the ratio 

< UU >N f =2 

is significantly smaller than one, one may expect a rather small value of Nf rtt \ In the 
functional integral, the dependence upon Nf arises from the fermion determinant part of 
the measure: setting all quark masses equal, one has 

dfi = dfi(A) (det(i Ip+m)) N f . (2) 

In other terms, it is a Dirac sea effect. In the quenched approximation, which is often 
used in lattice simulations, the fermion determinant is set equal to one and -R32 is exactly 
one. The same result also follows in the leading large N c expansion of QCD, since the 
determinant contributes to graphs with internal quark loops, which are subleading. Using 
chiral perturbation theory (CHPT), one can access the ratio 

5 < uu >( mu=md=ms= Q) 

-^32 = — — — ■ (A) 

< uu > (m u =m d =0,m s ^0) 

This ratio is different from R32, but it is also a measure of the influence of the fermion 
determinant in the evaluation of the quark condensate. Again, this ratio would be exactly 
one in the leading large N c expansion or in the quenched approximation, for any value of 
the strange quark mass m s . The point here is that the physical value of strange quark 
mass is sufficiently small compared to the scale of the chiral expansion A ~ 1 GeV, such 
that the chiral expansion in m s /A makes sense and, at the same time, m s is not that small 
so that -R32 will not trivially be close to one. In this paper, we will provide an estimate of 

^32- 

The plan of the paper is as follows. In sec. 2, the expression of -R32 in CHPT at 
order one-loop is given. This expression involves a single low-energy coupling constant, 
Le(fi), in the nomenclature of Gasser and LeutwylerQ. In that paper, Lg was simply 
assumed to be OZI suppressed. Here, we attempt a more careful estimate on the basis of 



a chiral sum rule. Analogous chiral sum rules were discussed in the recent literature [ 10 1 



and eventually provide very good precision [ 11 1. In sec. 3, a sum rule expression for Lq in 



terms of the correlation function ng(s) of two scalar currents uu + dd and ss is established 
and /mllg(s) is shown to satisfy a Weinberg- type sum rule. This will be an important 
constraint to our evaluation. The construction of the spectral function is discussed in 
sec. 4. Important ingredients are the pion and the kaon scalar form-factors which can be 
related to experimental information on pion-pion scattering using analyticity, unitarity, 
high-energy constraints as well as low-energy constraints from chiral symmetry. This was 



first performed in ref.[12]. Extension to the region of 1.5 GeV, where important resonance 



contribution is expected is then discussed. Finally, the result can be found in sec. 6. 

2. Ratio of quark condensates from CHPT 

Consider QCD in a limit where Nf quarks are exactly massless. We will consider the 
cases of Nf = 2 and Nf = 3 and assume, a priori, that chiral symmetry is spontaneously 
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broken in QCD in both cases and that the value of the condensate is sufficiently large 
also in both cases such that the conventional chiral expansion |l3|l ||9f applies. In nature, 
none of the quark masses m u , m s is exactly vanishing, but they are in an asymmetric 
configuration where m u , « m s (by a factor of twenty or so) and m s is itself sufficiently 
small compared to the scale of the chiral expansion A. Using this fact, one can express 
the ratio -R32 as an expansion in powers of m s . At chiral order 0(p 4 ), making use of the 
formulae of ref.||, one obtains 



R. 



32 



1 



m s B 
F 2 



32L 6 (/z) 
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16vH 



11 m s B 
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+ -ln- 
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At order 0(p 2 ) one has, 



m s B 



M 2 K+ + M 



Mi 



(4) 



(5) 



and this value can be used consistently in the equation above. The size of R32 depends on 
the value of a single low-energy coupling constant, Lq. 

The low-energy coupling constants of CHPT may be related to QCD correlation func- 
tions evaluated near zero momenta. This can be exploited, in particular for two-point 
functions, in order to express these constants in the form of sum rules using analyticity 
(and the fact that chiral correlators have non-singular short distance behaviour). A num- 



ber of these are exhibited in ref.|13|. A classic example concerns the coupling constant 
L10 which is related to the correlation function of two vector currents minus two axial 
currents. A very reasonable estimate for L10 can be obtained using simply the idea of 
vector meson dominance as well as Weinberg sum rules [ 14 1. Our aim is to estimate Lq 
along similar guidelines. 

One specific reason for interest in Lq is in connection with the Kaplan-Manohar 
transformation These authors observed that the effective lagrangian is left invari- 
ant under the following transformation of the quark mass matrix, 



M -> M + a 



/32So 
I F3 



(M*)~ l detM, 



(6) 



together with a transformation of certain low-energy constants. At chiral order 0(p 4 ) 
three coupling constants are affected, Lq, L-j and Lg, which get transformed as 



L-r — > 



a, 



L$ + 2a 



(7) 



One consequence is that using low-energy data alone, one can only determine combinations 
which are invariant under this transformation and not the individual values of Lq , L7 and 
Lg. These values are of some importance. In particular, the value of L$ determines the 
ratio of quark masses 2m s /(m M + md) beyond the leading chiral order p. It is therefore 
of interest to explore means of separately determining these constants (or at least one of 
them). 



3. Sum rule for L« 
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Consider the correlation function of the two scalar, isoscalar currents uu + dd and ss, 

n 6 (p 2 ) = iB 2 J d 4 xe tpx < T [(uu(x) + dd(x)) ss(0)] > c , (8) 

where the subscript c means that only connected graphs are to be retained. The fac- 
tor Bq 2 is introduced to simplify forthcoming expressions and furthermore makes 116 a 
renormalisation scale invariant object and Bq is defined as 

< uu > 



o 



lim ,,, 

mu=m d =m s =0 Fz. 



(9) 



For small momenta we can express Hq(p 2 ) using CHPT. In particular, at zero momentum, 
from CHPT at 0(p 4 ) one obtains, 



n 6 (0) = 64L 6 (/x)- 



1 



16vH 



„. (m s +m)B 4 (4m s + 2m)B 22 

2 m ^ — m ^ 

H 2 9 3/i 2 9. 



-0(m,m s ). (10) 



Here, and in the following, isospin breaking is neglected and we set m u = = m. This 
expression is at the basis of our sum rule estimate for Lq. The quark condensate ratio 
R32 has a very simple expression in terms of Hq. Combining equations (|j) and (|l0| ) one 
obtains 



R 



32 



1 



32vr 2 F2 



22 

i67r 2 n 6 (o) + — 



+ 0(m 



2 ) 

s)i 



(11) 



where barred quantities are to be taken in the limit m u = nid = 0. This relation (and 
eq.(|3)) can be recovered alternatively by noting that 



d 



dm. 



< uu + dd >= -B 2 U 6 (0), 



(12) 



and integrating this equation from m s = to its physical value using the CHPT expression 

It is possible to derive a lower bound on Lq and on Il^O) based on general properties 
of the QCD measure. At first, it is not very difficult to show that Hg(0) must be positive 
in the case of equal quark masses |j[6j]. Let Z be the partition function of euclidian QCD, 



dfi(A) e~ 5yM(A) det(z Jp A +M), M = diag(m, m, m s ) . 



We can express ^(0) in terms of Z as 

1 d 2 Z 



n 6 (o) 



1 dZ 1 dZ 



(13) 



(14) 



Z dmdm s Z dm Z dm s 
In the limit m s = m this can be written as an average of a manifestly positive quantity. 



lim n 6 (0) =< Tr 



1 

i If>+m 



< Tr 



1 

i Ip+m 



> >> 



(15) 
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where averages are defined as 

< >= i J dfi(A) e- SYM{A) det{i Ip A +M)0(A) . (16) 

This is because Tr (i Jj) +m) _1 can be shown to be real and the averaging is performed with 
respect to an integration measure which is real and positive in euclidian QCD (assuming 
the vacuum angle 9 = and a proper regularisation of the fermion determinant). We can 
now apply this result on the positivity of Ile(O) in conjunction with its one-loop expression, 
eq. (|To|) , setting all three quark masses there equal to the physical m s value. One gets 



2m s B 

log J - + 1 



+ 0{m s ) . (17) 



4608vr 2 

Ignoring higher loop corrections, this gives 

10 3 L 6 (Af„) > 0.35, 16vr 2 IT 6 (0) > 1.57 . (18) 
This shows, in particular, that the condensate must be a decreasing function of Nf. 

Another property, which will prove an important constraint for the sum rule estimate 
of Lq is that fle satisfies a Weinberg-type sum rule (WSR) 



o 



ImU 6 (s)ds = . (19) 



The proof is analogous to that of the ordinary Weinberg sum rules (see e.g. [fT^ ] ). The 
operators in the operator-product expansion at short distances must transform in the same 
way as (uu + dd) ss under the chiral group. The masses m u and being set equal to 
zero, the operator of lowest dimensionality having the correct transformation property is 
m s (uu + dd) and it has dimension four. Furthermore, a factor of (a s ) 2 is generated from 
the fact that all connected graphs contain at least two gluon lines. Taking into account 
the scale dependence of a s and that of -Bo in the perturbative region, we learn that IIq 
vanishes asymptotically faster than 1/q 2 , 

ft (n 2 \ - C<m s(uU + dd) > 1 

n6(9) "7[i^^ V' (20) 



with C a constant, which implies the WSR eq. (fL9|) . Furthermore, this behaviour ensures 
convergence of an unsubtracted dispersion relation for Ilg, so that one can express its value 
at zero as 

— , . If 00 , Imll 6 (s) , . 

n 6 o =-/ ds ^1 . 21 

vr Jo s 



While the WSR must be satisfied for arbitrary values of m s , it is clear from eq . fl20|) that 
convergence of the integral ( |l9| ) will be faster if m s = 0. In that situation, one expects the 
sum rule to be saturated in an energy interval of, say, [0 — 2] GeV, by analogy with the 
ordinary Weinberg sum rules. Experimental data is known for m s ^ 0, but while there 
may be large differences locally in imile(s) upon varying the value of m s (notably because 
of threshold effects), it is expected that these differences will be smoothed out to a large 
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extent in the integral. One therefore expects that the WSR will be also approximately 
saturated in a finite energy region, for the physical value of m s . 



Let us make some qualitative remarks on the practical significance of the sum rule. In 
the large N c counting, LTg is suppressed compared to a generic QCD correlation function, it 
is of order 0((N C )°) instead of 0(N C ). As a byproduct, one observes that the contribution 
of a single resonance is not enhanced by a factor of N c compared to the non-resonant 
background. In the large N c world, the coupling of a resonance to either uu + dd or ss 
will be suppressed. For a glueball, both couplings will be suppressed. The real world, 
in the scalar sector, seems to be quite different from these large N c considerations. For 
instance the /o(980) meson is found experimentally to be rather light, narrow, and it 
couples strongly to both KK and tttt channels in violation to the large N c expectation. 
As a consequence, one expects a strong contribution of the /o(980) to /mile. In order to 
satisfy the WSR ( |l~9| ) the contribution from the /o(980) has to be canceled by a higher 
energy contribution. It seems plausible that this will be resonance dominated as well. The 
particle-data book|l^| quotes several resonances in the 1.5 GeV region: the fo(1370), a 
rather wide resonance, the fo(1500) which is well defined and rather narrow and (possibly) 
the fo(1700). A first guess is that the sum rule (|l9| ) should be essentially satisfied from an 
interplay between the fo(980) and the fo(1500). 

The remaining problem is to estimate the couplings of these resonances to the scalar 
currents. This cannot be extracted directly from experiment because of the absence of 
a physical scalar iso-scalar source ( which is the same reason why Lq, L7, L% cannot be 
individually determined from low-energy experiments [|l9|l ). One way of getting around 
this difficulty, which was used in QCD sum rule estimates of the light quark masses [20] 



is to impose smooth matching of the resonance contribution with the low-energy domain, 
which is known from CHPT at leading order. This procedure can be checked to be a 
reasonable one in the case of vector currents. Implementation of this idea in the present 
context is discussed below. 

4. Construction of the spectral function 

A. The role of two-body channels 

In the construction of ImH^s) it is convenient to consider separately the two energy 
regions I) < s < 1 GeV 2 and II) s > 1 GeV 2 . Let us consider region I first. The only 
intermediate states allowed to contribute to Imlle(s) there are 2tt, 4tt and KK. When 
s << 1 GeV 2 , the 4tt contribution is suppressed by the chiral counting (being of order 
0(p 8 ) while the leading contribution is 0(p 4 )). Close to 1 GeV 2 , chiral counting is no 
longer effective, but it is found experimentally that the /o(980) has very little coupling 
to Air (in fact, no decay of the /o(980) into four pions has been observed yet []l8|| ). It is 
extremely likely, then, that the Air contribution to Imlle(s) is negligible in this whole 
energy range. As a result, the spectral function can be expressed in terms of the pion and 
of the kaon scalar form-factors. It is convenient to introduce the following normalisations, 

Fl{s) = T \fl < + Jd,7rV > Gl(s) = ^Vf < l Ia l 7r ° 7r0 > 
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F 2 (s) = — y/2 < 0\uu + dd\K + R- > G 2 {s) = —V2 < 0\ss\K + K~ > . (22) 
-Bo -So 

The values of these form-factors at s = are proportional to the derivatives of M% and 
Mjf with respect to the quark masses. At leading chiral order, one has 

Fi(0) = V6 Gi(0) = F 2 (0) = V2 G 2 (0) = V2. (23) 

One-loop corrections to these values can consistently be ignored because they are of the 
same order as the 0(p e ) contributions in eq, (|i~0|) . In the energy range I, the spectral 
function has the following expression 




AM 2 

ImMs) = F V \/ — -L F i (s)G*(s)e(s -AM?), (2 1) 

(where Mi = M n , M 2 = M K ). 

We consider now the energy region II. More approximations will have to be made 
in this region. We will work out the spectral function from several models in order to 
illustrate how the WSR can be satisfied. For the final purpose of evaluating the dispersive 



integral eq.(21) we will mainly rely on information from energy range I. As s increases a 
new two-body channel opens, rjrj. At some point, the Att channel will become important. 
Studies of tttt scattering suggest that this should happen at y^s ^1.4 GeV. This can be 
seen from the tttt inelasticity: below 1.4 GeV, inelasticity is found to be saturated to a 
good approximation by a single inelastic channel, KK (see e.g. fig. 7 of ref.|2l[]) and then 
one observes a strong onset of tttt — > Att. It is very likely that this is caused by the presence 
of the nearby scalar resonances /o(1370) and /q(1500) which were both observed to couple 



to four pion states [ 22 ] , 1 23 1 , [ 24 1 . Another experimental finding of these references is that 
the 47r system, in this energy region, likes to cluster into two resonances. This suggests 
that in an energy range sufficiently large to saturate the chiral sum rule, the contributions 
to the spectral function are either two-body channels or behave to a good approximation 
as quasi two-body channels. We will utilise below a model[|39| in which the An system is 
treated as an effective oo two-body channel. Correspondingly, we will introduce the scalar 
form- factors, 

F 3 (s) = — < 0\uu + dd\aa > G 3 (s) = — < 0\ss\aa > . (25) 
-Bo Bq 

This is certainly somewhat schematic as it is known that Air also clusters as an effective 
pp channel. In this model, furthermore, the r/rj channel is ignored. This is a questionable 
approximation, perhaps, although there is experimental evidence that the coupling of rjrj 
to tttt appears to be relatively suppressed p5j[ 26 1. In the quasi two-body approximation 



of multimeson channels one can express the spectral function in terms of form-factors 
Fi(s), . . . , F n (s) and Gi(s), . . . , G n (s) in the same way as eq.(p4|) except that the sum 
extends from 1 to n. Upon introducing effective two-body channels, one faces the difficulty 
that one can no longer rely on CHPT in order to determine the values of the form-factors 
at the origin. These values are needed in the construction to be described below. In 
practice, we will make the simple ansatz that these values are vanishing, 

F(0)=0, Gi(0) = 0, i>3. (26) 
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Evaluation of the scalar form- factors of the pion and the kaon was discussed in ref.[12|, 
based on a set of Muskhelishvili-Omnes equations. We review this evaluation and its 
extension below. 

B. Muskhelishvili-Omnes representation of scalar form-factors 

The form-factors Fi(s) which occur in the expression for the spectral function (p4|) 
(generalised to n channels) are themselves analytic functions everywhere in the complex 
plane except for a right-hand cut. Let Ty be the T-matrix elements which describe scat- 
tering among the various channels. A standard normalisation is adopted where the S and 
T matrices are related as 



4M? 

Sij = Sij + 2i<7 f 5 ry<7?0(s - 4M?)0(a - 4M'f), with a { = W - - 1 . (27) 

The discontinuity of the form-factors along the cut, generated from the two-body channels 
considered above, has the following form, 

n 

ImF^s) =Y J T* j {s)a j {s)F j {s)e{s - AM]) . (28) 

3=1 

We expect the form-factors to vanish asymptotically as [27] 

Fi(s) ~ l/s, s -» oo, (29) 
and therefore to satisfy an unsubtracted dispersion relation. Clearly, the approximation 



of quasi two-body channels cannot hold for arbitrarily large energies and eq.(28) is a 
reasonable approximation to the exact discontinuity only in a finite energy range. However, 
as we are interested in constructing Fi(s) in a finite energy region also, say below two GeV, 
the detailed behaviour of the spectral function at much higher energies is unimportant and 
we may as well assume that eq.(p8|) holds up to infinite energies, only requiring that the T- 
matrix behaves in a way that ensures the correct asymptotic decrease of the form- factors. 
Under these assumptions the form-factors must satisfy a set of coupled Muskhelishvili- 



Omnes [28]|29[](MO for short) singular integral equations, 



1 n r°° 1 Is' — 4M? 

F t {s) = - E / ds '~ T *A S ')\ 9(s' - UtJ)Fj(s') . (30) 

7T J ami s s V o J 

One observes that off-diagonal T-matrix elements are needed outside of the physical scat- 
tering region. Except in the one-channel case, this means that one not only needs physical 
scattering data but also a parametrisation model which allows for extrapolation. 

C. Asymptotic conditions on the T-matrix 

Let us now specify which asymptotic conditions are required from the T-matrix. Con- 
sider first the single channel case for which an analytic solution to the MO equation is 
available p8]] [ 29 1 , 



F(s) = P(s)Q(s), n(s) = exp 



, ds' - 1 7 S(s') 

7T J AMI \ S ~ S ) S 



(31) 
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where 5(s) is the scattering phase-shift and P(s) an arbitrary polynomial. Integrating by 
parts, it is simple to verify that, as s — > oo one has 

Ms) -» s~ l , I = -(6(00) - SOkMl)) . (32) 

7T 

Compatibility with the assumed high-energy behaviour of the form-factor is ensured pro- 
vided I > 1 . It is not difficult to see how this condition extends to the situation of n coupled 



channels [28], even though no analytical solution is known in general. Let us form a vector 
F of components (iq(s), . . . ,F n (s)); we learn from Muskhelishvili's book that there will 
be in general n independent solution vectors F a , a = 1, . . . n to the set of equations. Let 
us form an n x n matrix from these 

f(s) = (#(«),..., 400) • (33) 

All matrix elements of F(s) are analytic functions of s in the cut complex plane and the 
discontinuity across the cut can be formulated in matrix form, 

f(s + ie) = (l + 2»TE)F(s-ie), = 6ij(Ti(s)6(s - AM?) . (34) 

Taking the determinant of both sides, we obtain a one-dimensional discontinuity equation, 

f(s + ie) = D(s)f(s -ie), f = det F. (35) 

As f(s) is also an analytic function, this equation can be recast as a one-channel MO 
equation. As a consequence, the determinant of the solution matrix F can always be 
expressed in analytical form even though the individual entries are not known analytically. 
This is an interesting property which we have used as a check of the accuracy of our 
numerical calculations. It is easy to verify that, for a given value of the energy s with 
m < n channels being open, D(s) is the determinant of the m x m S-matrix so that it is 
a complex number of unit modulus, 

D(s) = exp(2iA(s)) . (36) 

Letting s go to infinity, det F behaves as an inverse power of s, so that the matrix F must 
be of the following form, 

lim f(s) = — C , (37) 

with C a constant n x n matrix with non- vanishing determinant. Taking the determinant 
of this equation and using the one-channel result one finds that 1/s asymptotic behaviour 
is ensured by the asymptotic condition, 

A(oo) - A(4M^) > nvr . (38) 

For instance, in the case of three channels, the T-matrix must be such that the sum of 
the three eigen-phase-shifts sum up to 37r (or more) when the energy goes to infinity. If 
the sum is exactly 37T, then the form-factors are uniquely determined at any energy from 
their values at zero. 
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D. Models of tttt scattering T- matrix 



In principle, tttt phase-shifts and inelasticities can be determined from di-pion produc- 
tion experiments, in which high-energy pions are scattered on proton targets (e.g. |pCf| ). 
A major source of information in this area so far, is from the high statistics experiment 



by the CERN-Munich collaboration [31] from which tttt S-matrix elements were extracted 
by a number of people[30|. Various determinations of S-wave phase-shifts are generally 
in reasonable mutual agreement below 1.4 GeV while marked differences are seen above. 
From these early analysis there was no clear evidence for resonances at 1.4 or 1.5 GeV. 
The CERN-Munich data themselves, however, are not incompatible with the presence of 



scalar resonances at these energies. This was demonstrated recently by Bugg et al. [ 32 1 
who obtained a good fit to the CERN-Munich data while constraining the S-matrix to 
have resonance poles and residues conforming to the PDG results. Unfortunately, it is 
not possible to use their parametrisation of the S-matrix for solving the MO equations 
because, on the one hand, it is not designed to satisfy the full set of two-body unitarity 
constraints and, moreover, the corresponding T-matrix parametrisation does not allow for 
extrapolation away from the physical scattering region. A set of tttt S-wave scattering 
phase-shifts and inelasticities was obtained recently |p3|, based on high statistics di-pion 
production data employing polarised proton targets pifl. In principle, polarisation infor- 
mation is extremely useful in reducing the problems of phase ambiguities. Two solutions 
consistent with unitarity were found, called up-flat and down- flat. We will consider the 
latter one only here, because on the one hand, it is in good agreement with earlier phase- 
shift determinations below 1.4 GeV and on the other hand, up- type solutions can usually 



be eliminated upon using the Roy equations [ 35 ] [3C] which encode crossing symmetry and 



high-energy constraints. This determination shows a marked resonance effect in the 1.4-1.5 
GeV region and thus appears as a good candidate for use in our sum rule analysis. One 
notes that, above 1.4 GeV, the 7T7t phase-shifts determined by [E3] and those determined 



in ref.[32| are not in good agreement. 



For the purpose of solving the MO equation system one further needs a T-matrix 
parametrisation allowing convenient (and reliable) extrapolation below physical thresh- 
olds. As pointed out in ref.[jl2| a useful check on the extrapolation of T\2 is to compare it 
to the chiral expansion in the region where the latter is valid. Close to s = one has, 

T 12 = -^s + 0(p*). (39) 

One-loop corrections to this result have been worked out [3^] [37]. A simple T-matrix model, 
which is very useful for performing checks of numerical calculations is that proposed by 
Truong and Willey[38[. This model has the property that the OM set of equations can be 
solved analytically. A somewhat more sophisticated T-matrix model, fitted to reproduce 



the tttt data of ref.j33[ was proposed in ref.[p9|. Fits with both 2-coupled and 3-coupled 
channels were performed. In this model, unitarity is ensured by solving a Lippman- 
Schwinger equation with a potential matrix chosen to have the following separable form, 

Vij( P , q) = J2 D 2 + »2 g 2 X lm • (40) 
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The T-matrix can be computed analytically and it can be checked to have the correct 
chiral magnitude at low-energy (in other terms it vanishes linearly with s and M%). It 
seems possible to adjust the parameters (and also the propagator) in order to reproduce 
exactly the correct T-matrix chiral expansion at 0(p 2 ) and even, we believe, at 0(p 4 ), but 
this has not yet been done. In this model, the OM equations must be solved numerically. 
For this purpose, we have developed an algorithm which is described in the appendix. 

The A and /i arrays in eq. (^) are constant parameters fitted to the data. In the case of 
3-coupled channels, the available data is not sufficiently constraining and several different 
sets of parameters can provide comparable fits. Two different sets of parameters were 
obtained in refs.|39| (called A and B) and two further sets in ref.(4^] (called E and F). Sets 
A, B and E generate fits with comparable x 2 with the set of data considered in ref . p9| . 
Set F has a good x 2 a t l° w energy only. Close to 1.4 GeV it has a very narrow resonance, 
which, perhaps, could be interpreted as a glueball. Although not producing a very good 
X 2 , the authors of ref.[^0| suggest that this scenario is not totally excluded by the data. 
The data which was used in these fits consists in a) The set of tttt phase shifts 5 n (E) and 
inelasticities ^(E) as determined in ref. |33|l in the energy range 0.6 < E n7T < 1.6 GeV and 
b)the set of phases 4>i2(E) of the tttt — > KK amplitudes from the particular experiment of 
Cohen et al.|41j|. One must keep in mind here that there is some discrepancy in the lower 



energy part between the result of this experiment and others, notably by Etkin et al. [ 42 1 
as far as the phase is concerned. This point is discussed in some detail by Au et al.pl|, 
whose K-matrix parametrisation could more easily reproduce the latter phase results. 
Regrettably, the absolute values of the tttt — > KK amplitudes, which are also available 



from experiment, were not included in the fits of refs.[39][40|. The various parameter sets 
differ in the behaviour of the phase-shift in the third channel, 63(E), which is unconstrained 
by experiment and also, to some extent, on the detailed structure of inelasticities. These 
differences, as we will see, will result in fairly different behaviour of the spectral functions 
as well, so that the sum rule (0) appears as an interesting theoretical constraint in this 
kind of analysis. 

The T-matrices generated from this model do not satisfy the asymptotic constraints 
eq. (|3^) neither in the 2-channel case nor (for any of the parameter sets discussed above) in 
the 3-channel case. Thus, they cannot be used up to infinite energies for our purposes. We 
must impose the proper asymptotic behaviour, i.e. that the eigen-phase shifts must sum 
to 2tt in the case of two channels and 3tt for three channels, by handf]. For this purpose, 
we have introduced a cutoff energy Eq. For E < Eq the T-matrix is computed from the 
model, while for E > Eq the phase-shifts are interpolated as follows 

5 W (E) = utt + (8«{Eq) - mr)f (JQ , 5i{E) = 5i(E )f ( , i > 2, (41) 

with n the number of channels and the cutoff function f(x) = 2/(1 + x m ). In practice, we 
have taken Eq = 1.5 GeV which insures a smooth raise of 5 n (E) and m = 3. Changing 
these parameters will modify the details of the shape of the spectral function in the higher 

3 These are the minimal asymptotic values which ensure existence of a solution. We will assume that a 
possible further raise above the minimal values can only occur for y^s > 2 GeV and will have no influence 
on lower energy results. 
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energy region. Inelasticities are computed from the model in the whole energy range. 
Other elements of the S-matrix can then be deduced from unitarity and continuity. 

6. Results 

A. Two-channel models 

We first calculate the scalar form-factors and the spectral function from irir — KK 
2-channel models. The result for the spectral function, using the T-matrix model of Au 
et al.0 |2lJ is shown in Fig.l, together with the result from the 2-channel version of the 
potential model of ref.|j{J]. Consider first the region y/s < 1 GeV: there, the spectral 
functions from the two models have the same sign and are rather similar in shape. For 
sum rule applications it is useful to introduce the following spectral function integrals in 
this energy region 

I n = 16vr / £±1 ds . 42 

Some numerical results for the integrals I n , n = 0, 1 are shown in table 1 below. The 
predictions from the two models are seen to differ by less than 10% for these quantities. 





Au 


KLL 


set A 


set E 


set F 




3.30 


3.09 


4.18 


3.38 


2.69 


h 


4.92 


4.73 


6.56 


5.74 


4.04 



Table 1: Values of the integrals Iq and I\ (see eq.(42)) from the 2-channel T-matrices 

The last three columns correspond to different parameter 



labeled as Au t 21] and KLL l 
sets of the 3-channel T-matrix model of Kaminski et al. 



We have also computed, for these two models, the low-energy observables associated 
with pion form-factors. In the neighbourhood of s = 0, one defines (we follow the notations 
of ref.10]) 

F 1 (s) = Fi{0) [l + ^<r 2 + c n s 2 + ..], (43) 
6 

and, similarly, for the matrix element of the ss current 



-Mg-Gi(s) = d F s[l + b A s + ...] 



(44) 



As shown in ref.|T^] the parameter dp is proportional to the derivative of F w with respect 
to the strange quark mass. Upon using the T-matrix from Au et al., we have verified 



that our calculation reproduces the results obtained previously [12] [43]. The numbers are 



4 As observed in ref. Jl2| upon using the values of the parameters at the precision given in the Au et al. 
paper, a spurious very narrow resonance appears close to the KK threshold: we removed this resonance 
by linearly interpolating the nn phase-shift on both sides. If not removed, the spectral function would be 
identical to that shown in the figure except at the very position of the resonance where a very narrow dip 
would appear. 
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displayed in table 2. We also show the results corresponding to the T-matrix from ref , p9| . 
The numbers quoted in the table correspond to an improved T-matrix where the irir 
phase-shift is constrained in the low-energy region ^/s < 0.6 GeV in order to match 
the predictions from CHPT at two-loops for the scattering length and the scattering 
range@]|5| i.e. ag = 0.21, 6g = 0.26M" " 



If we do not make this modification, the 



results from the T-matrix of ref. [39] would be only slightly different. For instance, one 
would have, < r 2 >J= 0.580 /to 2 , dp = 0.075 GeV~ 2 , reflecting the reasonable low-energy 
behaviour of the T-matrix in this particular model. These results are compatible with the 
known low-energy coupling constants from CHPT. At order one loop, the chiral expansion 
of < r 2 >J and dp involve the constants L4 and L5 



24 

< r * >"=J2 2L 4 ( M ) +J L 5 (/x) 



1 



32vr 2 



dF 



8M| 
F 2 



' M 2 1 M 2 
lo § -T + 7 lo § — T 

fl z 4 [A A 



L 4 (/i) 



256vr 2 



1 + log 



/x 2 



(45) 
(46) 



While L4 is not easily determined elsewhere, L5 can be extracted from the ratio of Fx/F n 
and this givesQ L 5 (M p ) = 1.4 ± 0.5 10~ 3 . Using eqs.(||) and the results from Table 2, 
one deduces 



L 4 (M P 



0.410" 



L 5 (M P 



110" 



(47) 



The value of L4 can be considered as a prediction, and the result for L5 appears to be 
compatible with Fx/F n . One must note, though, that this agreement might be somewhat 
fortuitous because the error of this determination is rather large. If we assume, for instance, 
10% relative errors on < r 2 >J and on dp, the resulting uncertainty on L5 would be 
AL 5 (M p ) = ±0.710~ 3 . 





< r 2 >l 


C-7T 


dF 


b A 


10 3 L 4 


10 3 L 5 


Au 


0.585 


10.50 


0.087 


3.29 


0.41 


0.81 


KLL 


0.605 


10.81 


0.076 


3.45 


0.36 


1.09 


set A 


0.609 


10.86 


0.134 


3.12 


0.61 


0.62 


set E 


0.583 


10.58 


0.144 


3.15 


0.66 


0.29 


set F 


0.653 


11.51 


0.047 


3.80 


0.23 


1.79 



Table 2: Results for some low-energy observables and corresponding chiral coupling con- 
stants at 0(p 4 ). < r 2 >JT is in fm~ 2 and other quantities in appropriate powers of GeV. 
The Erst two lines correspond to 2-channel T-matrices and the the last three lines to 
3-channel models. The labeling is the same as in Table 1. 

In the region E > 1 GeV now, the spectral functions from the two models differ consid- 
erably, see Fig.l. That corresponding to the model of ref.[^] exhibits a strong resonance 
effect. Its contribution to the integral goes in the sense of canceling the positive contribu- 
tion from the /o(980). This is qualitatively as expected from the WSR ( |l9| ) and indicates 
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that it seems possible to satisfy this constraint in a 2-channel model. Quantitatively, how- 
ever, if one uses the set of parameters of ref.|3j| without alteration, one finds that the 
contribution from the /o(1500) is somewhat too strong and overcompensates that of the 
/o(980). At any rate, in the 1.5 GeV region it is no longer a good approximation to retain 
only two channels in unitarity relations. Other two-body channels are open like rjrj, and 
experiment indicates significant coupling of the /o(1500) to the 4ir channel as well. We will 
now investigate how such additional channels can affect the results in an approximation 
of a single effective additional channel. 

B. Three-channel models 

We have computed the spectral function based on the 3-channel T-matrix model from 
ref.|3£| using several parameter sets determined in this reference and in a subsequent 
onepO[. Results are shown in Figs. 2-4. Fig. 2 corresponds to the parameter set A from 
ref.|^j (we recall that set B from this reference was discarded because T\2 has an unphys- 



ical low-energy pole in this case) and figs. 3,4 correspond to the sets E and F from ref.[40j 
respectively. Again, let us consider first the energy region y/s < 1 GeV: there, the spectral 
functions from the 3-channel models are comparable to those from the 2-channel ones. 
This can be seen, for instance, for the integrals Jo and I\ displayed in table 1. Also one 
can see from Table 2 that the result for < r 2 >J are rather stable. A somewhat less stable 
quantity is dp, the derivative of the pion form- factor of the ss current, which increases in 
the 3-channel model for parameter sets A and E. This results in a slight increase of L4 and 
a significant decrease of L5 which becomes too small and incompatible with Fx/F T . On 
the contrary, a large value of L5 emerges if one uses parameter set F. Keeping in mind the 
uncertainty in the determination of L5 from the scalar form- factors, none of the 3-channel 
parameter sets considered here is as satisfactory as the 2-channel models as far as the very 
low energy behaviour of the form- factors is concerned. 

Let us now consider the energy region E > 1 GeV. One observes from Figs. 2-4 that 
a variety of shapes can get generated from different parameter sets. Set F (see Fig. 4 ) 
has a negatively contributing resonance, but it is much too strong and does not obey the 
WSR eq. (lll|) . Set A has a positively contributing resonance and does not obey the WSR 
constraint either. Set E (see Fig.3) displays a more complicated structure: the /o(1500) 
has also a positive contribution but there is a wide negative contribution centered at 1.7 
GeV, and the WSR is approximately obeyed. 

C. Estimate of S 6 (0) 

The main conclusion from the above results is that while there seems to be reasonable 
agreement on the shape of the spectral function in the energy range yfs < 1 GeV, its 
structure above 1 GeV is subject to considerable uncertainty. In models with more than 
two coupled channels the parameters are not sufficiently constrained from the experimental 
7T7T and KK data. Another source of uncertainty in those models which, again, can be 
checked to affect the spectral function above one GeV concerns the values of the form- 
factors at the origin -Fj(O), Gj(0), i > 3 which are not given from chiral symmetry. At 
least, we have seen that there exist models which fit the data and can also accomodate 
the WSR constraint. 
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In order to calculate IT^O) from the spectral representation eq.(21) one needs, in 



principle, to know the spectral function both below and above 1 GeV. However, the lower 
energy range is expected to generate the largest contribution. Qualitative information on 
the high energy sector, such as the existence of the WSR constraint and the experimental 
position of the resonances is sufficient if one is not asking for a very high precision. Firstly, 
we expect the contribution to ^(O) from the range y/s > 1 GeV to be negative, giving 
the upper bound, 

n 6 (o) < j^h . (48) 

The most plausible scenario is that of a single resonance dominated contribution around 
yfs ~ 1.5 GeV. In this scenario, the following estimate of 5e(0) is valid, 

in which the WSR has been used. In this case, the correction from the higher energy range 
is approximately 30%. The other possibility is that several resonances, the /o(1500) and 
the /o(1700), are playing a role in the sum rules. If the two resonances make negative 
contributions one expects that the correction to 56 (0) will be smaller than 30% because 
the 1/s factor suppresses the /o(1700) contribution. If the contribution from the /o(1500) 
is positive the correction is even smaller (this was realised in one of the 3-channel models 
considered above). The last possibility is that of a negative /o(1500) and a positive 
/o(1700) in which case the contribution from the higher energy region will be largest, but 
simple estimates like ( fi"9| ) show that a 50% correction is a generous upper limit. This gives 
us a lower bound on 56(0), 

n 6 (o) > g^/i . (50) 

From these considerations and the numbers of table 1 we infer that the value of 56 (0) 
must lie in the following range 

2 < 16vr 2 5 6 (0) < 6 . (51) 



Using chiral perturbation theory to one loop, eq. (|1Q[) , this result can be recast into an 
estimate of the coupling constant Lq 

0.4 < W 3 L e (M ri ) < 0.8 . (52) 

We note that the bound ( |i"8| ) obtained in sec. 3 is satisfied. This number can be compared 
with the estimate of ref. M 

L 6 (M V ) = (0.0 ± 0.3) 10" 3 . (53) 

The central value there, is obtained from the assumption that the OZI rule applies. Indeed, 
the OZI rule implies that is identically equal to one, and inserting Lq(M ti ) = in eq.(||) 
one finds R32 = 0.96 which is very close to one. 

The value of Lq that we obtained from the sum rule implies the following result for 
the ratio of quark condensates -R32 

R 32 ~l- 0.54 ± 0.27 . (54) 
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In order to obtain this estimate, we have used for Lg the central value which emerges 
from the sum rule discussion above and, for the error, we have used the same value as that 
estimated in ref-S (see eq.(53) ). Within the substantial error band, the main observation 



is that the deviation of the quark condensate ratio from one is negative, and it seems to 
be rather large. 

7. Summary 

We started by noting that the sensitivity of the quark condensate on Nf can be tested 
by studying its variation as a function of the strange quark mass. This variation may be 
related to the correlation function IlQ^q 2 ). Two different expressions of 116(0) are equalled, 
one based on chiral perturbation theory and one which uses a dispersive representation. 
We discussed the spectral function which enters this dispersive integral. We first argued 
that ImllQ satisfies a Weinberg-type sum rule. This sum rule essentially relates resonance 
contributions from the two energy regions, < s < 1 GeV 2 and 1 < s < 4 GeV 2 . In the 
first energy region the spectral function can be expressed with good accuracy in terms of 
scalar form-factors of the pion and the kaon. In turn, these form-factors can be constructed 



from experimental scattering data on tttt and KK following the method of ref.[12|. The 
determination of the spectral function in the higher energy range is more uncertain. We 
considered the prediction from a model which treats the 4tt channel as an effective two- 
body channel. The influence of including this third channel in unitarity relations was 
found to have relatively minor influence on results in the lower energy region but has a 
strong influence on the region above one GeV. The dispersive integral, fortunately, receives 
its main contribution from the lower energy range. Using experimental information on the 
position of the resonances, as well as the WSR constraint, allowed us to obtain an estimate 

of n 6 (o). 

The conclusion of this analysis is that the properties of the /o(980) meson translate 
into a value of the coupling constant L§ which is significantly different from that expected 
from the OZI rule (or, alternatively, from large N c considerations). If one uses the central 
value obtained for Lq, one finds that the condensate uu decreases by as much as a factor 
of two as one decreases the mass of the strange quark mass from its physical value down 
to zero. A qualitatively similar behaviour is expected if one varies Nf from Nf = 2 to 
Nf = 3. This surprising result is possibly suggesting that Nf = 3 is not extremely far from 
a chiral phase transition point. One must bear in mind, however, that the relationship 
that we used between ng(O) and the condensate ratio receives corrections from two-loop 
CHPT. It remains to be seen whether these are significant or not. 
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Appendix: Numerical method 

The general idea for solving a linear integral equation is to approximate it by an 
ordinary linear system of equations by discretising the integral. The main difficulty in the 
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case of the MO equation is to handle the principal-value integral with high accuracy. Let 
us illustrate the method we have used on the one-channel MO equation, the generalisation 
to several channels is straightforward. First, one can transform the equation into one for 
the real part of the form-factor, R(s) = Re(F(s)), 



R(s) 



l 



f 



da' 



1 



vr J 4A/ 2 



-X(s')R(s') X(s') = taxxS(s'). 



(55) 



It is useful to split the integration region into several sub-intervals in order to accommodate 
fast variations of the integrand (we have used up to seven intervals in our numerical work) . 
Then, every sub-interval [a, b] is mapped to [—1,1] and the the quantity X(s')R(s') is 
expanded over a basis of Legendre polynomials. This allows us to perform the principal 
value integration in fl55|) using the exact formula p^]. 



-l u — z 



2Q L (z) 



(56) 



Here Ql(z) is the so-called Legendre function of the second kind. It is crucial, in order to 
ensure the success of the calculation, that it be computed to very high accuracy. An algo- 
rithm, based on using the recursion relations in the forward direction if \z\ < 1 and in the 
backward direction otherwise, proves adequate. One obtains a discretised approximation 
to the integral over [a, b] 



f 

J a 



N 



ds'^—X(s')R(s')^J2 W i 



1=1 



1 + 



2(s k - b) 



X{ Si )R{si) 



where 



-(a + b + (b 



a)ui), Wi[z] 



JV-l 

-Wi ^(2j + l)P.,(OQj 

j=0 



» 



(57) 



(58) 



and mat are the set of N Gauss-Legendre integration points (i.e. the zeros) of 

Pn(u) and w\, . . . ,wn are the associated set of weights. In the case where b = oo (last 
sub-interval), we use 

2a 

Si ~ 



and 



/oo 
ds'- 
a 



Sk 



-X(s')R(s') 



1 



2a 

Sk 



Ui 



N 



i 



2A 

Sk . 



X( Si )R( Si 

1 - Ui 



(59) 



(60) 



In this manner, the functional equation for the function R(s) gets transformed into a 
set of M linear equations for R(s\), . . . ,R(sm) where M = nN, n being the number of 
sub-intervals. We note that this is a homogeneous system which, strictly speaking, has no 
nontrivial solution unless the determinant vanishes. In practice, it does not exactly van- 
ish. It is only in the limit of iV — > oo, in fact, that the determinant vanishes. In addition, 
one wants to specify the value at zero R(0) and this generates one additional equation, 
which is non- homogeneous. A solution can be defined by dropping one of the homoge- 
neous equations. A numerically stable way of performing this, is to use the singular-value 
decomposition of the (M + 1) x M matrix of the linear equation system jj^]. 
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We have performed several checks of the numerical calculations: a) we have verified 



that upon using the T-matrix of the Truong-Willey type [38] the analytical result was 
accurately reproduced b)we have verified the stability of the result when varying the 
number of integration points up to several hundred points and c) we have also verified that 
the determinant of the n-independent solutions obtained numerically accurately reproduces 
the result which is known analytically (see sec. 4.C). 
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Figure 1: Spectral function computed from T-matrix models with two coupled channels. 
The solid line corresponds to the T-matrix of ref.^\], the dashed line corresponds to 
refMl 
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Figure 2: Spectral function computed from the 3-channel T-matrix model of ref. 3£ ]. The 
curve corresponds to the parameter set A of this reference 
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Figure 3: Spectral function computed from the 3-channel T-matrix model of ref. 4C ]. The 
curve corresponds to the parameter set E of this reference 
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